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We investigate theoretically the collective charge den- 
sity wave motion through an ensemble of small disordered 
Aharonov-Bohm rings. It is shown that the magnetic flux 
modulates the threshold field and the magnetoresistance with 
a half flux quantum periodicity $o/2 = h/2e, resulting from 
ensemble averaging over random scattering phases of multiple 
rings. The magnitude of the magnetoresistance oscillations 
decreases rapidly with increasing bias. This is consistent with 
recent experiments on NbSes in presence of columnar defects 
[Phys. Rev. Lett. 78, 919 (1997)]. 

PACS numbers: 03.65. B, 72.15.R, 72.15.N 



One of the first experimentally observed quantum phe- 
nomena in charge density wave (CDW) compounds is 
the oscillating magnetoresistance in NbSea in presence 
of small columnar defects, recently reported by Latyshev 
et al. 0. The collective response of the CDW to the 
Aharonov-Bohm flux trapped inside the columnar defects 
reflected a $o/2 = h/2e periodicity. The occurrence of 
this period was related to instantons ||^ in ring-shaped 
commensurate CDW conductors, where large scale quan- 
tum fluctuations of the CDW phase allow for macroscopic 
quantum tunneling between degenerate ground states. 
Related work ||] on CDW's in a ring geometry, however, 
predicts a periodicity of $0 due to the modulation of the 
amplitude of the order parameter. It is doubtful whether 
these models apply to an array of columnar defects in 
a planar film. Here we propose a different theoretical 
model, which is closer related to the actual experimental 
geometry, and can account for the observed effects. 

Figure [| shows schematically a planar film of one- 
dimensional CDW chains, containing a small hole 
threaded by magnetic flux. The CDW is characterized 
by a complex order parameter |A|exp(zx), where |A| is 
proportional to the amplitude of the density modulation 
and the phase x denotes its position. The size of the hole 
is smaller than both the longitudinal and transversal 
CDW coherence lengths. As recently shown by Arte- 
menko and Gleisberg Q , nonlinear screening of the phase 
distortions induced by a defect, leads to the formation of 
metallic islands surrounding the latter. Hence, in our 
model we include a small conducting region around the 
hole where the CDW order is destroyed. In the normal 
region, electrons can encircle the hole via different paths, 
as indicated by the dotted lines, each picking up a ran- 
dom scattering phase. Because of the strong anisotropy. 



the largest contribution will arise from the middle chains. 
Therefore, as a model calculation, we proceed with a one- 
dimensional treatment of the problem. 

Using microscopic equations for the sliding CDW mo- 
tion "over" a general scattering source, we show that 
the Aharonov-Bohm flux modulates the CDW threshold 
field. The periodicity $0/2 appears by ensemble averag- 
ing over random scattering phases, as is known from the 
Al'tshuler, Aronov, and Spivak (AAS)-theory 1^ for an 
ensemble of Aharonov-Bohm rings. Our results qualita- 
tively account for the amplitude of the magnetoresistance 
oscillations and its disappearance at higher bias. 

In the framework of the kinetic equations ^J, the 
motion of the quasi-particles and the condensate in a 
quasi one-dimensional CDW conductor can be described 
by the semiclassical Green functions g^^(x;i, i') where 
i = {R,A,K} and a,/3 = {1,2}. The retarded and 
advanced g'^ functions determine the excitation spec- 
trum and the Kcldysh component g^ describes the ki- 
netics of the system. The subscripts 1 and 2 refer to 
the right-, respectively, left-moving electrons at the two 
branches of the linearized energy spectrum. We incorpo- 
rate the effects of the Aharonov-Bohm ring in this for- 
malism, by imposing boundary conditions on the Green 
functions. 

The Green functions at the right (R) and left (L) of a 
small scattering source of characteristic length I < ^0 are 
related by 



g^;,^-(Mt)-ig^^Mt, 



(1) 



where M is the transfer matrix of the scatterer which 
satisfies the condition M^^ ~ cr^'M.eTz in order to en- 
sure current conservation. The transfer matrix M can 
generally be paramtrized as 



M 



(2) 



where rj^^ipo are scattering phases and A/(l + A) = R 
is the reflection probability of the scatterer j^. To in- 
vestigate the motion of the CDW in the vicinity of the 
scatterer we proceed in five steps. First we gauge away 
the phase x(a;, t) of the CDW by performing the unitary 
transformation 



(3) 



where U — exp^x'^'a- Then we apply a Fourier trans- 
formation with respect to the time difference in Wigner's 
representation 



1 



(4) 



with T — {t + t')/2. The third step is to restrict ourselves 
to the Keldysh Green function and integrate it over all 
energies. Using the identity / deg^de — Acr^ we obtain 
the boundary condition 



C J degj^ ~ J deg^C + 2h{xR -Xl)C = 0, (5) 

where C = U^^MtUi^, and xr.l = dtXR,L define the 
sliding CDW velocities at the right and left side of the 
scatterer, respectively. Next it is convenient to decom- 
pose the matrices into the unit matrix and the three Pauli 
matrices as g = gol + g ■ ^ and C = CqI + C ■ where 
9 = {91,92,93), ^ = {£i,'^2,C3) and a = (cri, 0-2, erg). 
After substitution into (0) we are left with the equations 



&(5ofl-5oL) + 2?i(Xi^-XL) = 



(6a) 



>Co J deig^ -g^)+^Cx J dt{f^ + 5l ) - 0. (6b) 

From the defenition of current 

eN{Q)vF 



I = 



J deTi-g'^ + eN{0)vphx, (7) 



where — e is the electron charge, vp is the Fermi velocity 
and A^(0) is the density of states at the Fermi energy, 
it is seen that the first equation just expresses the con- 
servation of current through the scatterer. In the quasi- 
stationary approximation, where we neglect the inertia 
of the CDW {dfx = 0) and the time dependence of the 
amplitude of the order parameter, we know from the self- 
consistency equation that 



J deg^ = 0> / deg^ 
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(8) 



where 7 is the dimensionless electron-phonon coupling 
constant ||. Substituting Eqs. (|) and (|) into Eq. (|^), 
we finally obtain 

XR-XL = 2770, (9a) 



^J■R- fJ-L - -;;-VRcos{x- (po) = 0. (9b) 

Here we have defined the electrochemical potential 4/i = 
— / deg^ , x = {xr + Xl)/2, and we have assumed equal 
amplitudes of the CDW order parameter at the right and 
left contacts |A|i^ = |A|l = |A|. The first equation rep- 
resents the correlation between the CDW's on the right 
and left hand side of the scatterer in equilibrium, and 
is necessary for a self-consistent treatment of pinning. 



The second equation is the desired equation of motion 
for the CDW. It is similar to the phenomenological sin- 
gle particle model for impurity pinning |^ with threshold 
potential fix given by fir — |A|a/R/27. Equation ( pb| ) 
allows us to calculate the conductance of a CDW moving 
over a small-size but arbitrary scattering source, which 
is characterized by its transfer matrix evaluated at the 
Fermi energy. 

Now that we have studied the problem of a general 
scatterer we proceed by calculating the transfer matrix 
of a disordered ring threaded by a Aharonov-Bohm mag- 
netic fiux. Following Biittiker the two junctions of 
the ring are described by the unitary scattering matrix 
S, which relates the outgoing to the incoming scattering 
amplitudes 



(10) 



where il± = — 2f3, and < (3 < i is a couphng pa- 

rameter describing the reflection at the entrance. We rep- 
resent the upper {+) and lower (— ) path around the ring 
with the transfer matrices N_|_ and N_, parametrized as 



Nh 



e'''± 
e-"'± 



(11) 



where ri± and (j)± are the scattering phases of the sepa- 
rate branches. For simplicity we consider here the clean 
limit without barriers in both arms, and associate only 
different scattering phases to different paths. A more 
general analysis will be presented elsewhere ||ll[. The 
phase shift of the electrons due to the magnetic flux is 
accounted for by the transformation N± exp(±iz9)N±, 
where ■& = 7r$/$o with ffux quantum <I>o = h/e. After 
some algebra we obtain the total transfer matrix M of 
the disordered ring 

AA/11 = Q'i cos 2(1) + n1 cos 2-d 

-4(1-/3) cos 2ri- Aif] sin 2r] 
AM22 = -(AAfii)* 
AM12 = - AAfai = Q'i cos 2(j) 



- cos 2-d + 4^1 - 2/3 cos 2ri, 



(12) 



where we have defined ij = {-q^ +77_)/2, = (77+ — 77_)/2 
and A = 4/3{exp(i77) cos(-!5-|-(/)) — exp(— ijy) cos(i? — 0)}. It 
is easily verified that this transfer matrix ensures current 



conservation. From Eq. (6b) we obtain Lp^ — 7r/2 and 
the self-consistency relation becomes 



Xfl - XL = 2 arctan ( 



tan T] 



n\ sin^ ^ + SIT? ( 



2/3 sin 277 



(13) 



which oscillates as a function of the flux and the scat- 
tering phases 77 and </>. The threshold field also depends 
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strongly on the magnetic flux and the scattering phases 
through 
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Mo, 



27 



(14) 



This expression is dominantly periodic in the flux quan- 
tum, but also contains higher harmonics from weak- 
localization paths. As expected for equal phases = 0, 
destructive interference at values $ = $o/2 mod $0 
enhances the total backscattering and thus the pinning 
force. For 7^ its behavior is more complex. 

So far we considered only a single Aharonov-Bohm 
ring. We now turn to the problem of an ensemble of 
rings in a CDW system. It is well known that the resis- 
tance of an ensemble of Aharonov-Bohm rings in series 
or parallel retains only the half flux quantum periodic- 
ity $0/2 This is due to the different scattering 
phases of subsequent rings. Interference effects of elec- 
trons traveling only half the ring circumference average 
out. However, if electrons encircle the rings just once, 
the phase difference of time- reversed paths is $o/2 in- 
dependent of the scattering properties of the individual 
rings. In our model the ensemble average comes about 
from averaging over multiple columnar defects. In the 
limit of strong pinning where \/R\A\/ {nhvpn) ^ 1, with 
large impurity potentials or a low impurity concentration 
n, it is known from the Fukuyama-Lee-Rice (FLR) model 
p3[ that the CDW adjust its phase to each defect in or- 
der to minimize the electrostatic Coulomb energy. As a 
consequence the net threshold field is proportional to the 
sum over all impurities. For uncorrelated defects in series 
the above formalism reproduces the FLR-model in this 
limit. If we assume a random distribution of scattering 
phases for the Aharanov-Bohm rings, we obtain the net 
threshold field Et by taking the ensemble average over 
the scattering phases 77 and </> of a single hole 



£;t($) = 



1 



dr] / d0^T($,?7,0), (15) 



where d is the average distance between the columnar de- 
fects. The result is shown in Fig. ^ for different values of 
the parameter /3. The average threshold field oscillates as 
a function of half the flux quantum $0/2 around a con- 
stant value. With increasing reflection /3 at the junctions 
the total threshold shifts upwards and higher harmonics 
become visible due to the increasing dwell-time of elec- 
trons in the rings. The appearance of a finite pinning 
force at zero field arising from the columnar defects is 
consistent with experiments. 

Above we have shown that the threshold field of an 
ensemble of columnar defects oscillates with a half fiux 
quantum period. To qualitatively understand the ef- 
fect of the flux dependent threshold held on the current- 
voltage characteristics, we distinguish between the high 
and low sliding velocity regimes. Near threshold the de- 
pinning and dynamics of CDW's can be described as a 



dynamical critical phenomenon |I 
IcDW obeys the scaling relation 

E — Ej" 



IcDW fx 



E^ 



The CDW current 



(16) 



where the critical exponent a takes the value 3/2 in 
mean- field ||l5| . Since the threshold field contains a flux- 
dependent term as in Eq. ([l5|), large oscillations in the 
magnetoresistance are expected with $o/2 periodicity. In 
the high velocity limit, however, the magnetoresistance 
oscillations disappear rapidly, since the CDW does not 
see the pinning potentials and their flux dependence. We 
remark that if we are allowed to neglect phase-slip pro- 
cesses, the potential drop over a single defect should be 
smaller than the gap |A|. The actual correlation between 
the phases left and right is probably smaller that in our 
model calculation, due to, for example, the finite size of 
the hole and the metallic region. If we take into account 
the energy dependence of the transfer matrix, it can be 
shown that the pinning force in Eq. (9b) is then modified 
by|A|^|A|exp(-//eo). 

In comparison with the experiments ||l| , this model ac- 
counts for the periodicity of the magnetoresistance oscil- 
lations and qualitatively for the decrease of the ampli- 
tude at large biases. The observed minimum at zero field 
may arise from spin-orbit interaction, which is known to 
determine the sign of the magnetoresistance oscillations 
in compounds with relatively large atomic- numbers 
This has not been taken into account presently and needs 
further investigation as well as an extension to multiple 
channels. Experimentally, measurements of the depen- 
dence of the threshold field on the magnetic flux should 
provide more insights into the validity of our model. 

We conclude by summarizing our results. We have 
derived microscopic equations for the non-linear sliding 
CDW motion over an arbitrary scatterer. It is shown 
that magnetic flux trapped in a columnar defect modu- 
lates the threshold field. In a low density ensemble of 
defects, averaging results in a half flux quantum period- 
icity of the magnetoresistance. The amplitude of oscilla- 
tion decreases rapidly with increasing bias. 
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FIG. 1. Schematic figure of a small columnar defect 
threaded by magnetic flux in a planar film of CDW chains. 
Electrons can encircle the hole through a metallic region, 
where the CDW order is assumed to be destroyed (shaded 
area). 

FIG. 2. Averaged threshold field Et in units of lA|/27ed 
as a function of magnetic flux for different values of 
/3 = 0.3; 0.35; 0.4; 0.45; 0.5 from top to bottom. 
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